Anomalous decay of quantum correlations of quantum dot qubits 
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We study the evolution of quantum correlations, quantified by the geometric discord, of two 
excitonic quantum dot qubits under the influence of the phonon environment. We show that the 
decay of these correlations differs substantially form the decay of entanglement. Instead of displaying 
sudden death type behavior, the geometric discord shows a tendency to undergo transitions between 
different types of decay, is sensitive to non-local phase factors, and may already be enhanced by 
weak environment-mediated interactions. Hence, two-qubit quantum correlations are more robust 
under decoherence processes, while showing a richer and more complex spectrum of behavior under 
unitary and non-unitary evolution. 
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The study of quantum correlations in realistic systems 
has, for a long time, been limited to the study of entan- 
glement, due to the fact that straightforward methods 
of calculating the amount of correlations in a two-qubit 
system have only been available for some entanglement 
measures, such as the concurrence 0, 0] or negativity 
0, 0] • Although entanglement itself is a very important 
resource for a number of applications [j| , including quan- 
tum computation, quantum cryptography, or teleporta- 
tion, separability (the lack of entanglement) does not au- 
tomatically exclude the presence of quantum correlations 
Q. This is in particular the reason why quantum com- 
putation models relying on mixed, separable (not entan- 
gled) states 0-0| are possible. 

The quantum discord 1(| 11 [ is a measure of quantum 
correlations (see, however Refs. [Ill E3 f° r Holevo-type 
and thermodynamic based measures) which captures cor- 
relations beyond entanglement; it is defined as the dif- 
ference of two classically equivalent formulas for mutual 
information and is non-negative. Due to the null volume 
of the set of zero-discord states [l3j |. discord measures 
are not expected to undergo sudden death which is char- 
acteristic for entanglement evolutions 14-16j]. The ge- 
ometric measure of the discord describes the amount of 
correlations in a quantum system by finding the minimal 
Hilbert- Schmidt distance to the set of zero-discord states 
17j . Recently, a lower [17| and an upper [18| bound 
on the geometric discord which can be calculated from 
a two-qubit density matrix have been found, which sub- 
stantially simplifies the problem of studying the evolution 
of the quantum discord and opens the path for a qualita- 
tive and quantitative description of the decay of quantum 
correlations in realistic open quantum systems. 

In this paper we study the evolution of the lower 
and upper bounds of the geometric discord of two ex- 
citon quantum dot (QD) qubits interacting with an open 
phonon environment in order to capture the physical 



aspects of decoherence effects on quantum correlations. 
The interactions present in the system and the resulting 
dynamics are well understood. The experimentally ob- 
served evolution on picosecond timescales [l9|, [2(| can be 
described by pure dephasing within the independent bo- 
son model |20l . l2l| . Super-Ohmic phonon spectral densi- 
ties 22 , 23[ (resulting from the actual form of the carrier- 
phonon coupling and the phonon density of states (23 |) 
are responsible for characteristic features of the dephas- 
ing which is non-exponential and always only partial. 
Furthermore, a finite distance between the QDs leads to 
a time-delayed interference of phonon wave packets trav- 
eling from the two QDs which induces an environment- 
mediated interaction between the dots (and a small en- 
hancement of the density matrix coherences) in addition 
to the exciton-exciton interaction present in the system. 
The fact that the complex evolution of this ensemble can 
be credibly described theoretically in combination with 
experimental accessibility to a wide range of pure ini- 
tial states (which are optically excited on femtosecond 
timescales) make this system ideal for the examination 
of the quantum-information properties of open quantum 
systems. 

The specific system under study consists of two QDs 
stacked on top of each other and interacting with a 
phonon reservoir. The single qubit states |0) and |1) 
correspond to an empty QD and an exciton excited in 
the dot, respectively. The system is described by the 
Hamiltonian 

H = ei(|lXl|®I)+ea(I® |1)(1|) + Ae(|l)(l| ® |1)(1|) 
+(|1)<1|®I)^/W(6t +6_ fe ) (1) 

+(i ® |ixi|) fk\b{ + b- h ) + ^b[b k , 

k k 

where I is the unit operator, are the transition en- 
ergies in the two subsystems, Ae is the biexcitonic shift 
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(1 2) 

due to the interaction between the subsystems, ft ' are 
system- reservoir coupling constants, bk^b\. are bosonic 
operators of the reservoir modes, and uJk are the corre- 
sponding energies (we put H= 1). 

Exciton wave functions are modeled by anisotropic 
Gaussians with the extension Zj_ in the xy plane and l z 
along z for the electron and hole in both dots. The cou- 
pling constants for the deformation potential coupling be- 
tween confined charges and longitudinal phonon modes 
have the form fju'^ = fk^ ±tk ~ d ^ 2 , where d is the distance 
between the dots and 



fk = 
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2qVc 



(cr c - cr h )e 



where V is the normalization volume of the bosonic reser- 
voir, k±,z are momentum components in the xy plane 
and along the z axis, cr Ci h are deformation potential con- 
stants for electrons and holes, c is the speed of longitudi- 
nal sound, and g is the crystal density. In our calculations 
we put <7 C — 8 eV, <7h = —1 cV, c = 5.1 nm/ps, g = 5360 
kg/m 3 (corresponding to GaAs), l± = 5 nm, and l z = 1 
nm. The distance between the dots is taken equal to 
d = 6 nm unless stated otherwise. 

The Hamiltonian ([1]) can be diagonalized exactly using 
the Weyl operator method [24], [25| and we find the evo- 
lution of the double QD subsystem following Ref. [26j . 
Since local unitary transformations do not change the 
amount of quantum correlations in the system, we can 
use the density matrix p(t) — e~ lHht p(t)e lHht , with 
H L = £i(|l)(l| ® I) + £ 2 (I® |1)(1|), where E z = e, - 
J2k \fk\ 2 /oJk instead of pit) in the study of the geomet- 
ric discord. Assuming a separable initial system-reservoir 
state, we find the evolution of the elements of the density 
matrix p(t), which are equal to 

W)]u = \M«\ m)] ij = \M^- iAti{t)+B,iit) , (2) 

with 

4)1 (*) = A m =Y^\9k? s\nu k t, (3a) 

4>a(t) = 4^| 3fc | 2 cos 2 (fc z d/2)sinw fc t- ASt, (3b) 

A 12 {t) = 0, (3c) 

A 13 {t) = A 2 3 = A 03 -A i, (3d) 

Boi(t) = B 2 = Bis = B 2 3 (3e) 

= ^|5fc| 2 (coscj fc t-l)(2?i fc + l), 

B 03 (t) = 4j]| fffc | 2 cos 2 (fc z d/2)(cos Wfc i-l) (3f) 
x(2n fe + l), 

Bja(t) = 4Soi-S 03 , (3g) 

where n k is the Bose distribution, gk — fk/^k, A£? = 
Ae-2Re2 fe w fe |5fe| 2 e ifc - d , and the indices i, j = 0, 1,2,3 
correspond to the two qubit states 1 00) , 1 01) , 1 10) , and 



1 11), respectively. For long times, the factors coswfct 
and sinwfc£ become quickly oscillating functions of k and 
their contribution averages to 0. Consequently, the phase 
damping factors decrease form their initial value of 
to a certain asymptotic value depending on the mate- 
rial parameters, system geometry and temperature, while 
the phase shift factors Aij affect the system evolution at 
small times and then average out to zero. As a result, the 
off-diagonal elements of the density matrix are reduced 
and the phase information is partly erased. 

Here, we are interested in the symmetric geometric dis- 
cord [l7( which may also be expressed in the form of a 
purity deficit (see Ref. [l8j |) 
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(Tv[p 2 AB ] - Tr[(M A ® M b )pab]) ■ 

(4) 

Specifically, the discord is formulated as a purity 
(quadratic Renyi entropy) deficit under global versus 
product local (Ma<S>Mb) von Neumann measurements. 
In case of two qubits, the lower bound on the discord is 
given by [I3| 



D' s = max(Tr[iT x ] - k x ,Tr[K y ] - k y ) . 



(5) 



where k x is the maximum eigenvalue of the matrix K x — 
\x){x\ + TT T and k y is the maximum eigenvalue of the 
matrix K y = \y)(y\+T T T. Here, \x) and \y) denote local 
Bloch vectors with components Xi = Tr[p j 4s( cr i ® I)] and 
yi = Tr[p^s(I (g> <7j)], and the elements of the correlation 
matrix T are given by T^j = Tr[pAB(&i ® <Jj)\ (stemming 
from the standard Bloch representation of a two-qubit 
density matrix pab)- The upper bound is given by [18| 
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Tr[^] - k v 



k x + Tr[L y ] - 
- Tr[L x ] - l x ) . 



(6) 



where l x and l y are the maximal eigenvalues of the ma- 
trices L x = \x){x\ + T\k y )(k y \T T ^ and L y = \y)(y\ + 
T T \k x )(k x \T, respectively, while \k x ) and \k y ) are the 
normalized eigenvectors corresponding to the eigenvalue 
k x of matrix K x and k y of matrix K y . In the case of 
symmetric two-qubit states, meaning pab = Pba, no 
minimization or maximization is needed in eqs. ((5|) and 

The upper and lower bounds often coincide, yielding 
the true value of the geometric discord. This is specifi- 
cally the case for pure states, Bell diagonal states, and 
states with vanishing local Bloch vectors, \x) = \y) = 
[Isj ]. Hence, it is straightforward to show that the geo- 
metric discord is equal to 1/2 for all maximally entangled 
two-qubit states [27j |. 

\il>) = v^|00) + V6e M |10) + v/fee^lOl) - ^e^+^ll)), 

(7) 

with 2a + 2b = 1. Furthermore, the discord of the Bell 
diagonal states (a = or b — 0) under phonon-induced 
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FIG. 1: Evolution of the X-state geometric discord at T — 77 
K for d = oo (a) and d — 6 nm (b); solid red line - |o— b\ = 0.3, 
dotted pink- \a — b\ = 0.15, and dashed blue - \a — b\ = 0. 



partial pure dephasing is equal to Ds{t) = 2\pij(t)\ 2 — 
1/2 exp[2Bj,j (t)], where ij = 12 for a = and 03 for 
6 = 0, and the appropriate forms of Bij(t) are given by 
Eqs. ( |3"g| ) and (|3f|), which, up to a normalization, yields 
the square of the concurrence. 
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FIG. 2: Evolution of geometric discord bounds at d — 6 nm 
for pure initial state (Q with a — l/v2 (a,b) at different tem- 
peratures for AE — for which the upper and lower bounds 
are equal (a) and at 77 K for AE — 6 ps _1 ; red solid line 
- lower bound, blue dashed - upper bound, pink dotted - 
AE — (b). Lower bound values corresponding to the three 
eigenvalues of the matrix K x (the minimum of which yields 
the geometric discord lower bound) at different temperatures 
for zero biexcitonic shift (c) and and at 77 K for AE — 6 ps _1 
(d). 

Let us first study the evolution of the mixed X-state, 
/ 



P = 



a ag 03 (t) \ 

b bg 12 (t) 

bg* 12 (t) b 

\ag* 03 (t) a J 



(8) 



which is significantly simpler (although hardly accessible 
experimentally), but already carries some of the prop- 
erties of the discord evolution of a pure initial state 



with all coherences present. The entanglement of such 
a state, measured by the concurrence, is equal to C{t) = 
max{0, 6|<7i2(t)| — a,a\g Q3 (t)\ — b} and is prone to sud- 
den death. The geometric discord is given by Ds{t) = 
(a\g 03 (t)\ - %i 2 (t)|) 2 + (a - b) 2 , if \a-b\< a\g 03 (t)\ + 
b\g 12 (t)\ and by D s (t) = 2a 2 \g m {t)\ 2 + 2fo 2 | 5l2 (i)| 2 , if 
\a-b\ > a|0o3(f)| + %12(*)| (for long times, if a ^ b). 
Hence, the discord will not undergo sudden-death-likc 
behavior, but, if a ^= b, it will display a transition be- 
tween two types of decay (there is no simple relation 
between the transition point and the point of entan- 
glement sudden death). The transition point coincides 
with the transition point between quantum and classi- 
cal decoherence indicated in Ref. 28]. This is illus- 
trated in Fig. rjj where the geometric discord of the state 
©, with gij(t) = exp(— iAy(t) + By(t)), is plotted as a 
function of time for different values of \a — b\. The left 
panel corresponds to infinitely distant dots, for which 
<?12 (t) — go3 (t) , and the transition which is induced by 
the smooth partial pure dephasing process is clearly vis- 
ible. On the right panel, a similar evolution of the dots 
separated by the distance d = 6 nm is shown, which 
additionally displays an enhancement of the geometric 
discord after a finite time. This effect is due to a posi- 
tive interference between phonon wave packets traveling 
from the two dots. Note, that the process is sufficient to 
induce quantum correlations in an initially uncorrelatcd 
state with a — b (which remains uncorrelatcd, if d — oo). 

The next step is to study the evolution of the lower 
and upper geometric discord bounds for an initial state 
([7|) with all non-zero coherences (a ^ and 6^0) under 
phonon-induced partial pure dephasing. For simplicity 
the studied state is taken with a = b = 1/4 (the local 
phases a and (3 do not change the values of the geomet- 
ric discord or either of its bounds). On Fig. [2] (a) the 
evolutions of the geometric discord are plotted at differ- 
ent temperatures for zero biexcitonic shift (the upper and 
lower bounds are equal in this case) . The 3 K curve shows 
a distinct point where the discord is not smooth, resem- 
bling the evolution of the X-state (|8]), which is absent at 
higher temperatures. To understand this, the evolutions 
of Di = ^[Kx] — ki, where ki are the three eigenvalues 
of the matrix K x (the minimum of Di yields the true 
lower bound of the geometric discord) for 3 and 77 K are 
plotted in Fig. [2] (c). At 3 K a crossing of two Di curves 
is observed which is caused by the positive interference 
of phonon wave packets, which is responsible for the en- 
hancement of the geometric discord for the X-state of eq. 
©. 

Fig. [5] (b) shows the evolution of the lower (red solid 
line) and upper (blue dashed line) bounds on the geo- 
metric discord for the same initial state at 77 K when 
the biexcitonic shift is nonzero. The biexcitonic shift in 
the absence of any decoherence processes causes a coher- 
ent oscillation between the initial, maximally entangled 
state, and the separable state |i/> SGp ) = 1 /2(|0) + |1)) ® 
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(|0> + |1)) (reached when AEt = (2n + 1)tt, where n 
is a natural number). Under phonon-induced pure de- 
phasing, the oscillations of entanglement are damped and 
display prolonged periods when the entanglement is zero 
(which is only possible when the damping process can 
lead to sudden death of entanglement), and are otherwise 
smooth while their amplitude is limited by the entangle- 
ment decay displayed by the zero-biexcitonic shift evolu- 
tion (2(| . The oscillations of the geometric discord, which 
without decoherence would mimic entanglement oscilla- 
tions, are substantially different. Firstly, the discord does 
not display sudden-death-type behavior, and approaches 
states with AEt = (2n + l)ir smoothly, reaching zero at 
t = (2n+ 1)tt/AE only, if inter-phonon-interference does 
not induce extra coherence in the system (at short times 
and/or long distances between the dots). This confirms 
the notion that, since the set of zero-discord (only clas- 
sically correlated) states has zero volume, decoherence 
processes will never lead to the sudden and permanent 
vanishing of the quantum discord fl3j ] . 

Furthermore, the evolution of the quantum discord 
induced by the biexcitonic shift leads to the situation, 
when the value of the geometric discord is greater than 
the corresponding zero-biexcitonic-shift value. This can 
be clearly seen in the inset of Fig. [2] (b), where both the 
lower and upper bound of the geometric discord exceed 
the zero-biexcitonic-shift curve (pink dotted line). This 
behavior is non-monotonous and symmetric (for constant 
decoherence) with respect to the maximally entangled 
points given by AEt = 2nn (for which the non-zero- 
biexcitonic-shift and zero-biexcitonic-shift lines have to 
coincide) . This shows that the dependence of the discord 
on quantum phase relations is non-trivial, and than non- 
local phase correlations may lead to an enhancement of 
quantum correlations in mixed states depending on the 
actual value of the phase factor. The fact that the lower 
and upper bounds on the geometric discord are differ- 
ent in this case is in agreement with predictions made in 
Ref. [I8[ . We surmise that the discontinuity of the upper 
bound and the sharp features of the lower bound of the 
discord are an artifact of the procedure of their genera- 
tion from the density matrix, while the actual curve of 
the geometric discord is continuous and smooth. Fig. [2] 
(d) shows the evolutions of the lower bound values, Di, 
corresponding to the three eigenvalues of the matrix K Xl 
the minimum of which yields the actual lower bound, to 
illustrate the origin of the irregular shape of the lower 
bound of the geometric discord. 

We have studied the evolution of the geometric dis- 
cord of a two QD qubit system under decoherence caused 
by the phonon environment, giving the lower and upper 
bounds on the discord where it was impossible to find its 
true value. We have shown that the discord does not dis- 
play sudden death type behavior, but reveals a number 
of characteristic features (which are not displayed by en- 
tanglement) under the influence of phonons, which cause 



a continuous and smooth partial pure dephasing process. 
Firstly, the evolution of the geometric discord often dis- 
plays a transition between different types of decay, which 
is particularly evident for initial entangled X-states, but 
has also been observed for maximally entangled pure 
states with all coherences present. The study of the evo- 
lution of the discord in these pure initial states showed 
the importance of non-local phase correlations; a shift 
in the phase can lead to the enhancement of quantum 
correlations in a mixed two-qubit state. Furthermore, 
the positive interference of phonon wave packets origi- 
nating from the two dots (interaction through a common 
reservoir) which is weak in the system and cannot gen- 
erate entanglement between separable states, does lead 
to the appearance of quantum correlations described by 
the discord. Hence, the study of the quantum discord 
in this realistic scenario shows, among other things, that 
quantum correlations are a common occurrence in mixed 
separable states. 
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